ScPCA Supplemental cell type report for SCPCL000040

This report contains a summary of cell type annotation results for library SCPCL000040. The goal of this report is to provide more detailed information about cell type annotation results as an initial evaluation of their quality and reliability.

Performing cell type annotation is an inherently challenging task with high levels of uncertainty, especially when using automated annotation methods. One way to address this is to use multiple cell type annotation approaches and compare the results, as we have begun to do in this report.

When multiple methods annotate a given cell as the same or similar cell type, this may qualitatively indicate a more robust annotation.

Note that the contents of this report will vary based on which cell type annotations are present. Further, be aware that different cell type annotation methods may assign different labels to the same or similar cell type (e.g., the different string representations B cell naive and Naive B cell), due to use of different underlying reference datasets. Please note that all cell type annotation reference datasets are derived from normal (not tumor) tissue.

This library contains the following cell type annotations:

• Annotations from SingleR, a reference-based approach (Looney et al. 2019). The BlueprintEncodeData dataset, obtained from the celldex package, was used for reference annotations.
• Annotations from CellAssign, a marker-gene-based approach (Zhang et al. 2019). Marker genes for cell types were obtained from PanglaoDB and compiled into a reference named brain-compartment. This reference includes the following organs and tissue compartments: brain, embryo, immune system, and vasculature.

Sample metadata

The below table summarizes clinical metadata for the sample associated with this library. Blue hyperlinks are present for any terms with an ontology term identifier associated with the displayed human readable value. These links will direct you to a web page with information about that ontology term identifier.

Sample ID	SCPCS000040
Diagnosis	Pilocytic astrocytoma
Subdiagnosis	NA
Tissue location	Thalamic
Disease timing	Initial diagnosis
Age	11
Sex	F
Organism	Homo sapiens
Sample type	patient tissue

Cell type Annotation Summary

The plots and tables included here detail the results from performing cell type annotation.

Statistics

SingleR cell type annotations

In this table, cells labeled “Unknown cell type” are those which SingleR pruned due to low-quality assignments. In the processed result files, these cells are labeled NA.

Annotated cell type	Number of cells	Percent of cells
CD4-positive, alpha-beta T cell	1551	26.4%
central memory CD8-positive, alpha-beta T cell	972	16.55%
common lymphoid progenitor	540	9.19%
central memory CD4-positive, alpha-beta T cell	536	9.12%
regulatory T cell	486	8.27%
effector memory CD8-positive, alpha-beta T cell	370	6.3%
dendritic cell	250	4.26%
CD8-positive, alpha-beta T cell	203	3.46%
plasma cell	161	2.74%
natural killer cell	157	2.67%
effector memory CD4-positive, alpha-beta T cell	156	2.66%
megakaryocyte-erythroid progenitor cell	120	2.04%
granulocyte monocyte progenitor cell	102	1.74%
hematopoietic stem cell	91	1.55%
hematopoietic multipotent progenitor cell	62	1.06%
erythrocyte	44	0.75%
neuron	29	0.49%
megakaryocyte	7	0.12%
class switched memory B cell	6	0.1%
keratinocyte	6	0.1%
melanocyte	5	0.09%
memory B cell	3	0.05%
naive B cell	3	0.05%
epithelial cell	2	0.03%
inflammatory macrophage	2	0.03%
microvascular endothelial cell	2	0.03%
monocyte	2	0.03%
common myeloid progenitor	1	0.02%
fat cell	1	0.02%
macrophage	1	0.02%
Unknown cell type	3	0.05%

CellAssign cell type annotations

In this table, cells labeled “Unknown cell type” are those which CellAssign could not confidently assign to a label in the reference list. In the processed result files, these cells are labeled "other".

Annotated cell type	Number of cells	Percent of cells
T helper cells	56	0.95%
T cells	42	0.72%
Nuocytes	12	0.2%
Microglia	9	0.15%
T regulatory cells	9	0.15%
Oligodendrocyte progenitor cells	6	0.1%
Monocytes	2	0.03%
Dendritic cells	1	0.02%
Endothelial cells (aorta)	1	0.02%
GABAergic neurons	1	0.02%
Gamma delta T cells	1	0.02%
Macrophages	1	0.02%
Myeloid-derived suppressor cells	1	0.02%
Natural killer T cells	1	0.02%
Neural stem/precursor cells	1	0.02%
NK cells	1	0.02%
Trophoblast cells	1	0.02%
Unknown cell type	5728	97.51%

UMAPs

In this section, we show UMAPs colored by clusters. Clusters were calculated using the graph-based Louvain algorithm with Jaccard weighting.

Next, we show UMAPs colored by cell types. For each cell typing method, we show a separate faceted UMAP. In each panel, cells that were assigned the given cell type label are colored, while all other cells are in grey.

For legibility, only the seven most common cell types are shown. All other cell types are grouped together and labeled “All remaining cell types” (not to be confused with “Unknown cell type” which represents cells that could not be classified).

Cell label comparison plots

This section displays heatmaps comparing cell labels from various methods.

We use the Jaccard similarity index to display the agreement between between pairs of labels assigned by different annotation methods.

The Jaccard index reflects the degree of overlap between the two labels and ranges from 0 to 1.

• If the labels are assigned to identical sets of cells, the Jaccard index will be 1.
• If the labels are assigned to completely non-overlapping sets of cells, the Jaccard index will be 0.

High agreement between methods qualitatively indicates higher confidence in the cell type annotation.

Unsupervised clustering

Here we show the labels from unsupervised clustering compared to cell type annotations. Cluster assignment was performed using the Louvain algorithm.

Automated annotations

This section displays a heatmap directly comparing SingleR and CellAssign cell type annotations. Note that due to different annotations references, these methods may use different names for similar cell types.

Quality assessments of automated annotations

SingleR annotations

To assess the quality of SingleR cell type annotations, we use the delta median statistic.

• Delta median is calculated for each cell as the difference between the SingleR score of the annotated cell type label and the median score of the other cell type labels in the reference dataset.
• Higher delta median values indicate higher quality cell type annotations.
  • Values can range from 0-1.
  • Note that there is no universal threshold for calling absolute high vs. low quality, as described in the SingleR book section on ‘Annotation diagnostics’.

You can interpret this plot as follows:

• Each point represents the delta median statistic of a given cell whose SingleR annotation is shown on the y-axis.
• The point style indicates SingleR’s quality assessment of the annotation:
  • High-quality cell annotations are shown as closed points.
  • Low-quality cell annotations are shown as open points. In other sections of this report, these cells are labeled as Unknown cell types.
  • For more information on how SingleR calculates annotation quality, please refer to this SingleR documentation.
• Diamonds represent the median of the delta median statistic specifically among high-quality annotations for the given cell type annotation.

CellAssign annotations

To assess the quality of CellAssign cell type annotations, we consider the probability associated with the annotated cell type. These probabilities are provided directly by CellAssign:

• CellAssign first calculates the probability of each cell being annotated as each cell type present in the reference.
• CellAssign then annotates cells by selecting the cell type with the highest probability among all cell types considered.
• These probabilities range from 0 to 1, with larger values indicating greater confidence in a given cell type label. We therefore expect reliable labels to have values close to 1.

The plot below shows the distribution of CellAssign-calculated probabilities for the final cell type labels. Line segments represent individual values that comprise each distribution.

For cell types with 2 or fewer labeled cells, only the individual value line segments are shown. Line segments are also taller for any cell type label with 5 or fewer cells.

Session Info

R session information

## ─ Session info ───────────────────────────────────────────────────────────────
##  setting  value
##  version  R version 4.4.1 (2024-06-14)
##  os       Ubuntu 22.04.4 LTS
##  system   x86_64, linux-gnu
##  ui       X11
##  language (EN)
##  collate  en_US.UTF-8
##  ctype    en_US.UTF-8
##  tz       Etc/UTC
##  date     2024-08-02
##  pandoc   2.9.2.1 @ /usr/bin/ (via rmarkdown)
## 
## ─ Packages ───────────────────────────────────────────────────────────────────
##  package              * version  date (UTC) lib source
##  abind                  1.4-5    2016-07-21 [1] RSPM (R 4.4.0)
##  beachmat               2.20.0   2024-04-30 [1] Bioconductor 3.19 (R 4.4.1)
##  beeswarm               0.4.0    2021-06-01 [1] RSPM (R 4.4.0)
##  Biobase              * 2.64.0   2024-04-30 [1] Bioconductor 3.19 (R 4.4.1)
##  BiocGenerics         * 0.50.0   2024-04-30 [1] Bioconductor 3.19 (R 4.4.1)
##  BiocNeighbors          1.22.0   2024-04-30 [1] Bioconductor 3.19 (R 4.4.1)
##  BiocParallel           1.38.0   2024-04-30 [1] Bioconductor 3.19 (R 4.4.1)
##  BiocSingular           1.20.0   2024-04-30 [1] Bioconductor 3.19 (R 4.4.1)
##  bluster                1.14.0   2024-04-30 [1] Bioconductor 3.19 (R 4.4.1)
##  bslib                  0.7.0    2024-03-29 [1] RSPM (R 4.4.0)
##  cachem                 1.1.0    2024-05-16 [1] RSPM (R 4.4.0)
##  Cairo                  1.6-2    2023-11-28 [1] RSPM (R 4.4.0)
##  circlize               0.4.16   2024-02-20 [1] RSPM (R 4.4.0)
##  cli                    3.6.3    2024-06-21 [1] RSPM (R 4.4.0)
##  clue                   0.3-65   2023-09-23 [1] RSPM (R 4.4.0)
##  cluster                2.1.6    2023-12-01 [2] CRAN (R 4.4.1)
##  codetools              0.2-20   2024-03-31 [2] CRAN (R 4.4.1)
##  colorspace             2.1-0    2023-01-23 [1] RSPM (R 4.4.0)
##  ComplexHeatmap         2.20.0   2024-04-30 [1] Bioconductor 3.19 (R 4.4.1)
##  crayon                 1.5.3    2024-06-20 [1] RSPM (R 4.4.0)
##  DelayedArray           0.30.1   2024-05-07 [1] Bioconductor 3.19 (R 4.4.1)
##  DelayedMatrixStats     1.26.0   2024-04-30 [1] Bioconductor 3.19 (R 4.4.1)
##  digest                 0.6.36   2024-06-23 [1] RSPM (R 4.4.0)
##  doParallel             1.0.17   2022-02-07 [1] RSPM (R 4.4.0)
##  dplyr                * 1.1.4    2023-11-17 [1] RSPM (R 4.4.0)
##  dqrng                  0.4.1    2024-05-28 [1] RSPM (R 4.4.0)
##  edgeR                  4.2.0    2024-04-30 [1] Bioconductor 3.19 (R 4.4.1)
##  evaluate               0.24.0   2024-06-10 [1] RSPM (R 4.4.0)
##  fansi                  1.0.6    2023-12-08 [1] RSPM (R 4.4.0)
##  farver                 2.1.2    2024-05-13 [1] RSPM (R 4.4.0)
##  fastmap                1.2.0    2024-05-15 [1] RSPM (R 4.4.0)
##  flexmix                2.3-19   2023-03-16 [1] RSPM (R 4.4.0)
##  forcats                1.0.0    2023-01-29 [1] RSPM (R 4.4.0)
##  foreach                1.5.2    2022-02-02 [1] RSPM (R 4.4.0)
##  generics               0.1.3    2022-07-05 [1] RSPM (R 4.4.0)
##  GenomeInfoDb         * 1.40.1   2024-05-24 [1] Bioconductor 3.19 (R 4.4.1)
##  GenomeInfoDbData       1.2.12   2024-06-25 [1] Bioconductor
##  GenomicRanges        * 1.56.1   2024-06-12 [1] Bioconductor 3.19 (R 4.4.1)
##  getopt                 1.20.4   2023-10-01 [1] RSPM (R 4.4.0)
##  GetoptLong             1.0.5    2020-12-15 [1] RSPM (R 4.4.0)
##  ggbeeswarm             0.7.2    2023-04-29 [1] RSPM (R 4.4.0)
##  ggforce                0.4.2    2024-02-19 [1] RSPM (R 4.4.0)
##  ggplot2              * 3.5.1    2024-04-23 [1] RSPM (R 4.4.0)
##  ggrepel                0.9.5    2024-01-10 [1] RSPM (R 4.4.0)
##  GlobalOptions          0.1.2    2020-06-10 [1] RSPM (R 4.4.0)
##  glue                   1.7.0    2024-01-09 [1] RSPM (R 4.4.0)
##  gridExtra              2.3      2017-09-09 [1] RSPM (R 4.4.0)
##  gtable                 0.3.5    2024-04-22 [1] RSPM (R 4.4.0)
##  highr                  0.11     2024-05-26 [1] RSPM (R 4.4.0)
##  hms                    1.1.3    2023-03-21 [1] RSPM (R 4.4.0)
##  htmltools              0.5.8.1  2024-04-04 [1] RSPM (R 4.4.0)
##  httr                   1.4.7    2023-08-15 [1] RSPM (R 4.4.0)
##  igraph                 2.0.3    2024-03-13 [1] RSPM (R 4.4.0)
##  IRanges              * 2.38.1   2024-07-03 [1] Bioconductor 3.19 (R 4.4.1)
##  irlba                  2.3.5.1  2022-10-03 [1] RSPM (R 4.4.0)
##  iterators              1.0.14   2022-02-05 [1] RSPM (R 4.4.0)
##  jquerylib              0.1.4    2021-04-26 [1] RSPM (R 4.4.0)
##  jsonlite               1.8.8    2023-12-04 [1] RSPM (R 4.4.0)
##  kableExtra             1.4.0    2024-01-24 [1] RSPM (R 4.4.0)
##  knitr                  1.48     2024-07-07 [1] RSPM (R 4.4.0)
##  labeling               0.4.3    2023-08-29 [1] RSPM (R 4.4.0)
##  lattice                0.22-6   2024-03-20 [2] CRAN (R 4.4.1)
##  lifecycle              1.0.4    2023-11-07 [1] RSPM (R 4.4.0)
##  limma                  3.60.3   2024-06-16 [1] Bioconductor 3.19 (R 4.4.1)
##  locfit                 1.5-9.10 2024-06-24 [1] RSPM (R 4.4.0)
##  lubridate              1.9.3    2023-09-27 [1] RSPM (R 4.4.0)
##  magrittr               2.0.3    2022-03-30 [1] RSPM (R 4.4.0)
##  MASS                   7.3-60.2 2024-04-26 [1] RSPM (R 4.4.1)
##  Matrix                 1.7-0    2024-04-26 [2] CRAN (R 4.4.1)
##  MatrixGenerics       * 1.16.0   2024-04-30 [1] Bioconductor 3.19 (R 4.4.1)
##  matrixStats          * 1.3.0    2024-04-11 [1] RSPM (R 4.4.0)
##  metapod                1.12.0   2024-04-30 [1] Bioconductor 3.19 (R 4.4.1)
##  miQC                   1.12.0   2024-04-30 [1] Bioconductor 3.19 (R 4.4.1)
##  modeltools             0.2-23   2020-03-05 [1] RSPM (R 4.4.0)
##  munsell                0.5.1    2024-04-01 [1] RSPM (R 4.4.0)
##  nnet                   7.3-19   2023-05-03 [2] CRAN (R 4.4.1)
##  optparse             * 1.7.5    2024-04-16 [1] RSPM (R 4.4.0)
##  pillar                 1.9.0    2023-03-22 [1] RSPM (R 4.4.0)
##  pkgconfig              2.0.3    2019-09-22 [1] RSPM (R 4.4.0)
##  png                    0.1-8    2022-11-29 [1] RSPM (R 4.4.0)
##  polyclip               1.10-6   2023-09-27 [1] RSPM (R 4.4.0)
##  purrr                  1.0.2    2023-08-10 [1] RSPM (R 4.4.0)
##  R6                     2.5.1    2021-08-19 [1] RSPM (R 4.4.0)
##  RColorBrewer           1.1-3    2022-04-03 [1] RSPM (R 4.4.0)
##  Rcpp                   1.0.12   2024-01-09 [1] RSPM (R 4.4.0)
##  readr                  2.1.5    2024-01-10 [1] RSPM (R 4.4.0)
##  rjson                  0.2.21   2022-01-09 [1] RSPM (R 4.4.0)
##  rlang                  1.1.4    2024-06-04 [1] RSPM (R 4.4.0)
##  rmarkdown              2.27     2024-05-17 [1] RSPM (R 4.4.0)
##  rstudioapi             0.16.0   2024-03-24 [1] RSPM (R 4.4.0)
##  rsvd                   1.0.5    2021-04-16 [1] RSPM (R 4.4.0)
##  S4Arrays               1.4.1    2024-05-20 [1] Bioconductor 3.19 (R 4.4.1)
##  S4Vectors            * 0.42.1   2024-07-03 [1] Bioconductor 3.19 (R 4.4.1)
##  sass                   0.4.9    2024-03-15 [1] RSPM (R 4.4.0)
##  ScaledMatrix           1.12.0   2024-04-30 [1] Bioconductor 3.19 (R 4.4.1)
##  scales                 1.3.0    2023-11-28 [1] RSPM (R 4.4.0)
##  scater                 1.32.0   2024-04-30 [1] Bioconductor 3.19 (R 4.4.1)
##  scpcaTools             0.4.0    2024-07-08 [1] local
##  scran                  1.32.0   2024-04-30 [1] Bioconductor 3.19 (R 4.4.1)
##  scuttle                1.14.0   2024-04-30 [1] Bioconductor 3.19 (R 4.4.1)
##  sessioninfo            1.2.2    2021-12-06 [1] RSPM (R 4.4.0)
##  shape                  1.4.6.1  2024-02-23 [1] RSPM (R 4.4.0)
##  SingleCellExperiment * 1.26.0   2024-04-30 [1] Bioconductor 3.19 (R 4.4.1)
##  SparseArray            1.4.8    2024-05-24 [1] Bioconductor 3.19 (R 4.4.1)
##  sparseMatrixStats      1.16.0   2024-04-30 [1] Bioconductor 3.19 (R 4.4.1)
##  statmod                1.5.0    2023-01-06 [1] RSPM (R 4.4.0)
##  stringi                1.8.4    2024-05-06 [1] RSPM (R 4.4.0)
##  stringr                1.5.1    2023-11-14 [1] RSPM (R 4.4.0)
##  SummarizedExperiment * 1.34.0   2024-05-01 [1] Bioconductor 3.19 (R 4.4.1)
##  svglite                2.1.3    2023-12-08 [1] RSPM (R 4.4.0)
##  systemfonts            1.1.0    2024-05-15 [1] RSPM (R 4.4.0)
##  tibble                 3.2.1    2023-03-20 [1] RSPM (R 4.4.0)
##  tidyr                  1.3.1    2024-01-24 [1] RSPM (R 4.4.0)
##  tidyselect             1.2.1    2024-03-11 [1] RSPM (R 4.4.0)
##  timechange             0.3.0    2024-01-18 [1] RSPM (R 4.4.0)
##  tweenr                 2.0.3    2024-02-26 [1] RSPM (R 4.4.0)
##  tzdb                   0.4.0    2023-05-12 [1] RSPM (R 4.4.0)
##  UCSC.utils             1.0.0    2024-04-30 [1] Bioconductor 3.19 (R 4.4.1)
##  utf8                   1.2.4    2023-10-22 [1] RSPM (R 4.4.0)
##  vctrs                  0.6.5    2023-12-01 [1] RSPM (R 4.4.0)
##  vipor                  0.4.7    2023-12-18 [1] RSPM (R 4.4.0)
##  viridis                0.6.5    2024-01-29 [1] RSPM (R 4.4.0)
##  viridisLite            0.4.2    2023-05-02 [1] RSPM (R 4.4.0)
##  withr                  3.0.0    2024-01-16 [1] RSPM (R 4.4.0)
##  xfun                   0.45     2024-06-16 [1] RSPM (R 4.4.0)
##  xml2                   1.3.6    2023-12-04 [1] RSPM (R 4.4.0)
##  XVector                0.44.0   2024-04-30 [1] Bioconductor 3.19 (R 4.4.1)
##  yaml                   2.3.9    2024-07-05 [1] RSPM (R 4.4.0)
##  zlibbioc               1.50.0   2024-04-30 [1] Bioconductor 3.19 (R 4.4.1)
## 
##  [1] /usr/local/lib/R/site-library
##  [2] /usr/local/lib/R/library
## 
## ──────────────────────────────────────────────────────────────────────────────